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ANALYTICAL STUDY OF 
RIGIDIZED FIBROUS MATERIALS 


BY 

B. Walter Rosen, Debal Bagchi 
And John J. Kib ler 


1. INTRODUCTION 

Thermal protection during the re-entry of the space shuttle 
vehicle into the earth's atmosphere requires the development of an 
external surface insulation material having low density and 
moderately high temperature capability. This material must satis- 
fy two important requirements . It must be able to provide thermal 
insulation and it must be strong enough to prevent failure due to 
thermal, mechanical, acoustric and vibratory loads. The material 
must also offer satisfactory resistance to damage and degradation 
from ground handling environments such as fuel spillage, humidity 
and salt spray. To satisfy these broad requirements, insulation 
systems have been developed having fibers bonded together with a 
small amount of binder material. The materials have a very high 
void content and a very low density. The materials may be described 
as rigidized fibrous materials having a random orientation of fibers 
in three dimensional space. Methods for analyzing such materials 
have not been available. 

The present study involved the development of analytical material 
models suitable for the evaluation of the thermal and mechanical 
characteristics of various rigidized fibrous materials. The model 
treats the material as a space frame structure and includes effects 
due to fiber curvature and joint fixity. The present material model 
can be used to assess the relative effect of the material variables 
upon the strength under combined states of mechanical and thermal 
loadings. Thus, the relative merits of different materials (exist- 
ing or potential) can be studied. In the present task, attention 
was confined to mechanical and thermal loading. The effect of 
acoustic, vibratory or shock loads on the structural integrity was 
beyond the scope of the present research activity. 


2. MECHANICS OF FIBROUS MATERIALS 

The literature on the mechanics of fibrous materials is limited 
and largely related to the mechanics of paper. The material propert- 
ies of materials having a fibrous structure were treated by Cox [1] . 
He studied the effect of orientation of the fibers on stiffness and 
strength. For the planar case of fiber orientation distribution, it 
was shown that all possible combinations of elastic moduli may be 
represented by arrays of four sets of parallel fibers in appropriate 
ratios. Cox compared the results of the analysis with experimental 
results- for resin bonded fibrous filled materials and reported re- 
asonable agreement. 



For materials like paper or non-woven textiles, the micro- 
scopic distribution of fibers shows an irregular geometrical pattern 
Corte [2] and Dodson [3] have treated the strength and stiffness of 
this class of materials considering the statistical geometry of ran- 
dom fiber networks. The structure of this type of material is essen 
tially tv/o-dimensional ; i.e., no fiber traverses the thickness of 
the sheet from top to bottom. Thus, the geometrical problem of de- 
scribing the structure of this network is one in two dimensions. The 
basic approach to a quantitative description of a network of random 
lines representing the geometry of paper was outlined by Kallmes 
and Corte [4]. They expressed certain material characteristics in 
terms of measureable quantities including the fiber length, fiber 
width, number of fibers per unit area, eccentricity of orientation, 
etc. The dependent statistical measures of .the material included 
number of fiber crossings in the sheet, mean number of fiber cross- 
ings per fiber and the mean free fiber length. These quantities 
influence the strength and the elastic behavior of the network. 

Dodson [5] presented a statistical theory of the deformation 
of an inhomogeneous material under simple conditions of stress. 

The results of the analysis give statistical predictions of the ob- 
servable dilatations and the overall modulus in terms of the statis- 
tical properties of the material. Good agreement with experimental 
observation is reported for porous materials in the form of flat, 
fibrous networks. 

Rosen and Bagchi [6] studied behavior of a three dimensional 
fiber array by. using a space frame material model. That model con- 
sists of trusses, arranged in orthogonal planes, and having elements 
which differ in various aspects of their geometry including fiber 
diameter, nodal point spacing, fiber orientation, etc. Each indi- 
vidual space frame element can be analyzed when subject to applied 
loads to determine maximum stresses in the space frame members. 

The failure analysis then provided a relationship between maximum 
allowable applied stress and the various geometric parameters of 
the specified element. 


3. METHOD OF ANALYSIS 


The Reusable Surface Insulation (RSI) materials studied are 
rigidized fiber networks. The analysis of RSI materials is based 
upon the modeling of the actual material by a space frame struc- 
ture in which each section of fiber between nodes_ is treated as a 
beam having-some initial curvature. The model utilizes a series 



of different, regular arrays of trusses, having various orienta- 
tions, to represent the actual three dimensional fiber. array. 

This may be regarded as representing a continuous statistical dis- 
tribution function by a series of discrete delta functions of dif- 
fering magnitude (i.e., approximating a curve with a bar graph). 
Methods of translating statistical distributions of fibers into 
regular arrays have received some attention in the study of paper 
(Ref. 1 - 6) . The analysis of RSI materials requires treatment of 
nonhomogeneous materials. When this material is subject to any 
loading, the average deformation of the material depends both on 
the physical and geometrical variabilities of its constituents. 

The theory of deformation of heterogeneous media in the case of 
multiphase bodies considering physical variability from phase to 
phase has been studied by a number of authors (Ref. 7 - 11). The 
work that is associated with the class of materials termed as paper 
was discussed in the previous section and is concerned basically 
with the effect of geometrical variability of the constituent fibers 
upon the elastic properties of the material. 

The present work follows the concept that the three-dimensional 
random network of fibers can be represented by a basic set of truss 
elements oriented in different planes, as introduced by Rosen and 
Bagchi [6], Results obtained- from the analysis of that model show 
encouraging agreement when compared with experimental values. This 
has motivated the present development of a general model to evalu- 
ate the characteristic of RSI material exhibiting both geometrical 
and material variability. The geometrical and material variabilities 
of the actual material are represented by a number of trusses oriented 
in different planes with variable material and geometrical properties. 
The concept of parallel trusses is introduced to represent such geo- 
metric factors as orientation, fiber aspect ratio, eccentricity, 
volume fraction, etc. 

Basically, the approach is to represent the continuous 
statistical distribution function defining fiber geometry by a 
series of elements having discrete values of each of the vari- 
ables. This may be viewed physically as an approximate model, 
or it may be viewed mathematically as a discrete representation 
of a continuous function. 

The required analysis consists of three major aspects. First, 
the governing equations for the material model must be defined and 
a method of solution established. Second, a procedure for relat- 
ing the model parameters to the actual material geometry must be 
defined. Third, a failure criteria must be incorporated into the 
stress analysis procedure. Details are discussed below. 


3.1 Basic Assumptions 

The space model consists of an assembly of different trusses, 
groups of which are oriented parallel to each of five different 
planes (Figure 2) . Any of these five planes can have any number 
of trusses whose interaction with one another stems from the 
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constraint that all have the same average strains. The analysis 
of this model when subject to a combined state of stress is 
based on the following assumptions: 

1) Each of the trusses lying in any plane has the same 
state of average strains (which means that their edges 
have the same displacements) . 

2) The fibers within a single truss are of the same length, 
diameter, eccentricity and the same inclinations - 
oriented symmetrically with respect to X., axis (see 
Figure 1) . The material variability is accounted for by 
having parallel trusses with different values of para- 
meters assigned to different trusses. 


3.2 Stiffness Formulation 

Based on the above assumptions , the strain component in the 
fiber direction of any truss is given by, 

= B.e (1) 

where e - material strain components ; 

e f - fiber strains; 

B - transformation matrix. 

Similarly, the stress in the fiber of any truss is given 

by. 


o f = E.B.e (2) 

where - fiber stresses; 

E - effective fiber stiffness 

(including eccentricity effects) . 

At any stage of loading, the increase in the strain energy 
of each fiber link, 5V, of the model due to virtual displacements, 
6e, can be found: 

6V = f (e, E, B, G) .6 e (3) 

where G defines fiber geometry. The total strain energy is ob- 
tained by summation- over all links in a truss. Similarly, the 
additional external work done due to virtual displacement can 
be written as: 


6W = a . <5e 


( 4 ) 


where a are the applied stresses. 



Equating the total increase in strain energy to the addi- 
tional external work to obtain the equilibrium condition, we 
can write 


o = C* . e (5) 

where C* are the effective elastic properties of the material 
at any strain level. Details of the derivation are given in the 
Appendix A. 

The above equations define the nonlinear stress-strain 
relations as a function of the space frame geometry and fiber 
stresses. A computer program has been developed and checked 
out to compute these stress-strain curves for a relatively 
arbitrary space frame geometry. 


3.3 Elastic Symmetry 

If the internal composition of a material possesses sym- 
metry of any kind, then symmetry can be observed in its elastic 
properties. This will occur because the elastic properties 
are identical to the directions of symmetry which are developed 
in the body. The connection between symmetry of construction 
and elastic symmetry can be formulated in the following way: a 
material, in regard to its physical properties, has the same 
kind of symmetry as its crystallographic form. This principle 
can be expanded to include bodies which are not crystalline, 
but which possess a symmetry of structure. However, elastic 
symmetry is usually more extensive than geometric symmetry; in 
addition to the equivalent directions which coincide with the 
symmetric directions of the structure, other directions exist 
for which elastic properties are identical. 

If an anisotropic body possesses an elastic symmetry, 
then the equations of the generalized Hooke's law are simpli- 
fied. The effective elastic constants of the material model 
under study shows three planes of elastic symmetry. The mate- 
rial is orthotropic and has nine elastic constants. 


3.4 Thermal Stress 

For fibrous insulation materials , the principal mode of 
heat transfer at elevated temperature is by thermal radiation, 
see [13]. However, due to the random orientation of the fiber, 
the path of radiation will be obstructed causing a temperature 
gradient within the material. Since the material is under a 
uniform state of strain, the temperature gradient will cause 
thermal stresses locally within the material. These strain 



components may be treated as initial strains; the stress-strain 
law given by equation (5) can be written as: 


a = C*. (e - e T ) 


( 6 ) 


3.5 Determination of the Geometry of the Truss 

One of the problems of the space model is to obtain the 
regular geometrical orientation of truss models to represent 
the random oriented texture of the actual material. Corte [2] 
and Dodson [3] , based on certain statistical data for the in- 
dependent variables , provided the shape and area of the regular 
polygons formed by the truss members in two dimensions. In 
the present analysis, the model is defined by associating each 
fiber with one of the five reference planes. Each of these 
five groups of fibers is then regarded as a planar array and 
then statistical methods can be applied to define the geometry 
in each plane. As an example, the geometric variables may be 
regarded as independent quantities in each plane and they can 
be measured experimentally from enlarged photographs for that 
plane to determine a single variable frequency distribution 
function. Depending on the shape of the frequency curve, and 
the coefficient of variation, the area under the curve can be 
divided into an equal number of intervals to represent the 
continuous function with a fixed number of discrete variables. 
Detailed study of representing the random geometry with a 
fixed number of discrete variables is included in Appendix C. 


3.6 Failure Criteria 

When the model ''is subject to a combined state of stress, 
failure of the random oriented fibrous composite material may 
be viewed physically either as an accumulation of internal 
fractures caused by the breaking of fibers in any plane, or 
as a result of debonding of sufficient number of joints with 
the increase in applied stress. The class of rigidized fibrous 
materials dealt with here does not show any debonding effect 
during the experiment. As such, our main objective is to develop 
a failure criteria based on the fracture of individual fibers. 

The statistical models proposed in the study of fracture 
or failure mechanics take as a starting point Griffith's theory 
[14] which states that the reason for the difference between 
the calculated strengths of materials and the actual observed 
values resides in the fact that there exist flaws in the body 
which will weaken it. This means that there will be a distri- 
bution of strengths in a given specimen in the sense that a 
different amount of force will be needed to fracture the speci- 
men at one or another point. If one assumes that the flaws 
are distributed at random with a certain density per unit volume, 
the strength of a given specimen is determined by the weakest 



point in the specimen. This is the basis of statistical theory 
of brittle fracture of fibrous materials and other solids on 
which much theoretical work has been done [15]. However, the 
problem of distribution of strength of bundles of threads 
(fibers) [16] cannot be treated as a problem in the distribu- 
tion of smallest values because the complete bundle will not 
fail when its weakest member is broken. 

In a real problem, the situation is probably neither that 
of elements in series (weakest link) nor in parallel (bundle 
theory) , but elements distributed in rather complicated arrange- 
ment. Recently, the statistical theory of brittle fracture has 
been extended [17, 18] to define the failure of a fibrous body 
as the failure of the weakest link of a chain where each of 
the links represents a bundle of fibers connected in parallel. 

None of these theories described above are sufficient to 
describe the physical phenomenon of failure of RSI material 
with a random oriented fibrous structure. The problem is more 
complicated because of the fact that the stress field is non- 
uniform and the distribution of fibers is random. To get an 
exact solution, it is necessary to derive a joint density - 
distribution function describing the variability of geometrical 
and physical parameters and then correlate the nonuniform stress 
field from the above distribution function with the strength 
distribution. First of all, it is difficult to derive a joint 
distribution describing all variabilities. Secondly, after 
the initial of crack in some region, it will be necessary to 
modify the distribution function at each state of loading. An 
exact treatment of this problem based on existing knowledge of 
the statistical theory of brittle fracture poses an insurmount- 
able problem. • \ 

To simplify this problem, it is postulated in the present 
analysis that the material model consisting of an assemblage 
of truss elements is adequate to describe approximately the 
deformation behavior of the actual material in the pre and post 
crack region. With this postulation and further assumption 
that the truss will fail when any of its members will reach 
the maximum allowable stress, the probability of the failure 
of the complete model [19] based on the weakest element theory 
at a stress equal to or less than a° is given by 


n 

P(cr°) =1-71 [1 - F (S . ) ] ( 7) 

i=l 

Where F(S^) is the cumulative probability of failure of i th 
truss at a stress equal to or less than S-^ and n is the total 
number of trusses. 

Using the Weibull's distribution, the value of F(S^) 
is given by: 



where the values of a and 3 depend on the properties of RSI 
material . 

However, this formula is not satisfactory to describe 
failure conditions in the subcritical region after the initia- 
tion of a crack due to the purtubation of the stress field in 
the uncracked members of the trusses. To overcome this diffi- 
culty, a numerical approach is adopted based on the following 
definition of failure criteria to study deformation in the sub- 
critical region. As the loading is increased on a model, at 
some stage, fracture is initiated in a plane by the breaking 
of some fibers in a particular truss. Loading will be redistri- 
buted as a result of this break into the remaining unbroken 
fibers in different trusses. This mav lead to failure of more 
fibers with andther redistribution of stresses. However, at 
some point a stable configuration may be attained, arresting 
the growth of damage as a result of the redistribution. During 
this stable condition, loading can be increased, until eventually 
an unstable situation will arise when any further growth of 
fracture will cause a complete failure of the truss model due 
to redistribution. 

Therefore , failure of the present model can be viewed as 
shown in Figure (3 ) as a result of initiation of growth of frac- 
ture and finally due to the occurrence of an unstable situation 
when the reduced stiffness of the structure cannot absorb the 
energy released due to fracture. 

It should be mentioned here that depending on the stress 
condition, failure of any truss may be either due to the break 
of any fiber under tension and compression or due to local 
buckling of any fiber under compression. Both modes of failure 
can be detected using this numerical approach. 


4. PARAMETRIC EVALUATION OF ANALYSIS MODEL 

Parametric analyses have been performed using the space 
frame analysis program to obtain a feel for the sensitivity of 
the results to the many input parameters. Among the variables 
considered are : the effect of fiber joint fixity, the effect 

of distribution of both fiber angle and strength on tensile 
strength, compressive strength, shear strength, combined shear 
and tension, and thermal stresses. 

Figures 4 and 5 were picked at random to demonstrate the 
nature of the analysis. Figures 4 and 5 show the results for 
the RSI model with fiber edges free to rotate at the joints. 

Figure 4 gives the stress-strain plot of the model in some ar- 
bitrary units with one truss in each plane, while Figure 5 shows 
similar results for a model with five trusses in each plane, in- 
corporating the variability of Young's modulus and ultimate stress 
of the fiber in the trusses within each plane. 



Appendix B presents a detailed discussion of the results 
of the other cases considered during this study. The results 
clearly demonstrate that the stress geometry and fiber properties 
significantly effect the predicted RSI strength. 


5. RSI MATERIAL EVALUATION 

Before modeling RSI material as a space frame, photomicro- 
graphs of RSI were studied to determine the properties of the 
spatial-distribution of fibers, fiber diameter, initial curva- 
ture and fiber length to diameter ratio. Two types of RSI mate- 
rial were studied, the differences being primarily the diameters 
of the fibers. The material composed of smaller diameter fibers 
will hereafter be designated RSI-SDF while the material composed 
of larger diameter fibers will be designated RSI-LDF. Typical 
photomicrographs of these materials were supplied to MSC by 
NASA-LRC, two of which are reproduced in Figure 1. Given a least 
squares fit of a distribution function to the experimentally 
measured values, one then determined a discrete set of values to 
represent the continuous distribution function. Details of how 
the above was accomplished, along with the necessary analytical 
expressions are given in Appendix C. The remainder of this 
section will center on the results obtained from the photomicro- 
graphs of RSI and the corresponding analytical models which were 
developed. 

The angular distribution functions used to represent the 
spatial distribution of fibers were 


f(0) = ~ (1 + e cos 2 9) 

(9) 

f(0) 1 + e cos 2 ~ 

it + 2 e sin 29 

(10) 


Where equation (9) represents a simple elliptical distribution 
function and equation (10) provides symmetry in each quadrant, 
reflecting how one might expect a non-oriented fiber distribution 
to shift under processing pressures. 

The coordinate axes for An RSI tile are shown in Figure 6. 
The fiber distribution data and least squares best fit distribu- 
tion function using equation (9) are given in Figures 7 through 
10. As can be seen, a considerable amount of scatter existed in 
the data and equation (9) does not give a wholly satisfying fit 
to the data for the RSI-LDF, X-Z plane (Figure 7) . Note that 
for each group of points only one angle was measured. The spread 
of points about each angle was merely a convenient way of repre- 
senting multiple readings. 

The first set of Figures, 7 through 10, indicate that some 
eccentricity in the fiber distribution existed for the through- 
the-thickness planes, X-Z and Y-Z , for the RSI-LDF material, but 
that the in-plane photographs indicated a uniform, i.e. , circu- 
lar, distribution for both materials. 



The photomicrographs of the RSI-SDF material showed mainly 
fiber ends present in the through-the-thickness plane. Additional 
photomicrographs of this material were obtained and confirm the 
initial photographs. Due to the lack of distinct in-plane fiber 
definition in the through-the-thickness planes, the fiber distri- 
bution for RSI-SDF was not obtained in these directions. 

Figures 11 through 13 give the results of fitting the fiber 
distribution data to equation (10) . The best fit to the RSI-SDF 
material in the X-Y plane results again in a uniform distribu- 
tion, i.e. , e = 0.0, and is not included with the figures. Note 
that the best fit curves using equation (10) appear to provide 
a better fit to the data than that given by equation (9) . From 
the above data one recognizes that values of 6 can range from 
0 to 25 degress along with values of e varying from 0 to over 
0.2. With this information, limits can be set for the sensiti- 
vity studies such that the study results will fall in the realm 
of "real world" materials. 

In addition to fiber orientation data, the photomicro- 
graphs were used to determine fiber diameters of the various 
materials. Roughly one hundred measurements were taken from 
photomicrographs of the various planes of material. The results 
of these measurements are shown in Figures 14 and 15. It is 
clear from these figures that the RSI-SDF material fiber dia- 
meters were essentially uniformly distributed between 1 and 3 
micrometers in diameter, while the RSI-LDF material fiber dia- 
meters were roughly uniformly distributed between 2 and 6 mi- 
crometers in diameter. The RSI-LDF material, however, possessed 
several fibers in each plane which were up to 10-12 micrometers 
in diameter, which gave the photomicrograph a rather nonuniform 
diameter variation appearance. The RSI-SDF material, on the 
other hand, had little variation in diameter above 3 micro- 
meters, giving the photomicrographs a very uniform appearance. 

Fiber length to diameter measurements were virtually im- 
possible to make from the photographs since very few fibers 
possessed more than one intersection point in the photomicro- 
graphs, due to masking by other fibers, and due to the fibers 
simply extending out of view of the photomicrograph. At the 
other extreme, some fibers had intersections visible every 
one to two fiber diameters . 


5.1 RSI Model Definitions 

The analytical models used to represent RSI material are 
described in this section. The results of fitting the experi- 
mental fiber distribution data to elliptical distributions have 
been used to define the range of variables to be considered, 



since one cannot expect data from relatively small number of 
photomicrographs to completely characterize the real material. 

In an effort to fully cover the range and to determine the 
sensitivity of the results of the variables, mean fiber angles 
up to +40 degrees from vertical and eccentricities up to 0.30 
have been considered. 

Using the techniques described in Appendix C, the proper 
choices of angular positions for the fiber pairs in each plane 
have been made. As an example, Figures 16 and 17 give the dis- 
tribution function, f(0), and locations of fiber angles for various 
numbers of fibers in a plane for a mean fiber angle of zero degrees 
and eccentricity of zero (uniform distribution) and 0.100. The 
number of fibers per plane were varied to determine the number of 
fiber pairs required to give an adequate representation of the 
material. The results of this study will be discussed in the 
following section. 

Figure 18 presents the effect of varying the preferred fiber 
orientation for equation (10) from zero to forty degrees. Also 
shown are the discrete fiber locations for each preferred orienta- 
tion, assuming eight fiber pairs in a plane and an eccentricity of 
0.30. While the differences in fiber position are slight, some 
rather pronounced effects on the strength and modulus of the mate- 
rial were found, as described in the following section. 

Four planes were used in the models, located at positions of 
0, +45 and 90 degrees in the 1-2 plane. The differences in results 
obtained from the photomicrographs did not substantiate a signifi- 
cant difference in distribution between the xz and yz planes, hence 
the fiber distributions were assumed the same in all planes. 

Fiber length to diameter ratio, 1 / d, and eccentricity, e f , 

were held constant for each model, but were varied to determine 
the sensitivity of the results to each variable. Eccentricity 
was varied from zero (straight fibers) to 3.0 and 5,/d was varied 
from 5 to 15. The fiber volume fraction of all models was held 
constant at 10%. 

5 2 

A fiber modulus of elasticity of 1.378 x 10 MN/m 

(20 x 10 6 PSI) was assumed for all cases, and for the purpose 
of determining geometric effects, the fibers were all assumed 

to have a strength of 344.5MN/m^ (50KSI). In all cases ten- 
sile loading was applied to the 33 direction, or parallel 
to the planes containing the fibers. 


5.2 Discussion of the Results of the Geometry Sensitivity Study 

5.2.1 Number of Fiber Pairs/Plane Required 

The first problem which was investigated was to de- 
termine the number of fiber pairs per plane required to give 
reasonable results in the space frame analysis program. Since 



the cost per analysis was a strong function of the number of fiber 
pairs per plane, it was desirable to minimize this quantity, but 
not sacrifice confidence in the results. To provide an answer to 
this problem, a uniform distribution of fibers was assumed, with 
an eccentricity of 1.0. Values of SL/d were varied from 5 to 15 
and cases were run with 5, 8 and 10 fiber pairs per plane. The 
results are plotted in Figure 19. As can be seen, the number 
of fiber pairs per plane did not affect the predicted modulus of 
elasticity, but did alter the strain and stress at first fiber 
break . 

It should be pointed out that one pair of fibers at +9 
form a truss, which is assigned to a given plane. There can 
be any number of trusses assigned to a plane, and five planes 
in the model. Thus 8 fiber pairs per plane indicates that the 
model had 8 trusses defined (16 fibers) in each plane. 

The results appear to be fairly constant for the 8 and 
10 fiber pair cases, indicating that additional fiber pairs would 
not significantly alter the results. It- was concluded that using 
8 fiber pairs per plane would give sufficient accuracy for the 
purpose of this study. Consequently, eight fiber pairs per plane 
were used in the remaining analyses. For cases where multiple 
H/d or fiber eccentricities would be included in one model, the 
above conclusion would be altered since additional variables would 
be introduced in a single model, thus more trusses per plane 
could be required. 

It should be emphasized at this point that the choice of 
a uniform fiber strength for each model resulted in the first 
fiber break initiating complete failureof the model. Stable con- 
figurations after the first fiber break were only obtained when 
fiber strengths were variable in the model. The results of the 
analyses have been presented as : the elastic modulus in the direc- 
tion of loading, E 33 , measured by the stress and strain at first 
fiber break; the ultimate tensile stress, UTS 33 , which was simply 
the stress at first fiber break; and the strain at first fiber 
break; £y 33 , which was the strain at which instability of the model 
was initiated. 


5.2.2 Effect of Fiber Diameter 

The analyses formulated for the space frame analysis 
program indicate that, at least for zero fiber curvature, the re- 
sults are independent of fiber diameter and are a function of 
the fiber length to diameter ratio. In addition, since the model 



was based on the analysis of a representative unit volume, when 
one specified fiber angle and length to diameter ratio, the pro- 
gram calculates the appropriate fiber diameter for each fiber, 
consistent with the dimensions of the unit volume. It is possible 
to generate a model having the same fiber diameters and %/d. for 
all fibers by specifying the length of truss for each fiber pair, 
provided the lengths remain within the unit volume. 

As a check on the program and to determine if different 
fiber diameters affected the results, cases were run of the same 
model with constant fiber diameters of 5 micrometers, 0.318 mm, 
and variable fiber diameters. A fiber distribution with 9=0 and 
eccentricity equal to 0.2, and fiber eccentricity (curvature de- 
fined by initial displacement at the center of each fiber in terms 
of fiber diameters) equal to 1.00, keeping the results consistent 
for a fiber volume fraction of 0.10. The results are plotted in 
Figure 20. As can be seen, the results are indeed independent 
of fiber diameter and the length to diameter ratio is the pri- 
mary variable. 


5.2.3 Effects of Fiber Distribution Eccentricity 

With the above preliminary considerations satisfied, 
the more important effects of the distribution of fibers in the 
material were considered. The effects of eccentricity of fiber 
distribution, orientation of the fiber distribution, eccentri- 
city of the individual fibers and fiber length to diameter ratio 
were investigated. 

First, the eccentricity of the fiber distribution was 
varied, keeping the mean angle equal to zero and the eccentricity 
of the individual fibers equal to 1.0. The results are plotted 
in Figures 21 and 22 for various values of %/d. Note that 
eccentricity of the fiber distribution has virtually no effect 
on the strain at first fiber break although modulus of elasticity 
and applied stress at first fiber break increase significantly with 
increased distribution eccentricity. The results also show a con- 
sistent, although not large, effect of fiber length to diameter 
ratio. One would initially expect lower values of %/d to yield 
higher strengths since the truss members are becoming shorter for 
a given diameter, and the results indicate the opposite effect, 
Figure 22. However, the fiber eccentricity was held constant for 
these cases, which means that all center span initial deflections 
were one fiber diameter. This implies that for the smaller %/d 
values the fibers had much higher initial curvature than for the 
higher %/d values, and the higher initial curvature resulted in 
decreased strength, hence, the results are consistent. 


5.2.4 Effects of Orientation of Fiber Distribution 
On Properties 

The effects of the orientation of the fiber distribution 
defined by equation (10) on modulus of elasticity, strength, and 
strain at first fiber break are shown in Figures 23 through 25 , 
respectively for various values of eccentricity of the fiber 
distribution, holding %/d equal to 10 and eccentricity of the 
individual fiber equal to 1.0. 



Note that as the eccentricity of the distribution in- 
creases, the effect of orientation of the distribution becomes 
significant on the resulting strength and modulus of elasticity. 

As much as a 20% change in modulus and strength occurs for an 
eccentricity of 0.3, varing orientation from 0 to 40 degrees. 

Note, however, that the orientation of the fiber distribution 
has essentially no effect on strain to first fiber break, Figure 
25, regardless of the value of the eccentricity of the distribu- 
tion (i.e., a maximum strain change of only 0.6% for an eccen- 
tricity of 0.3). 

Also plotted on Figures 23 and 24 are the distribu- 
tion eccentricities and orientations obtained from fitting 
equation (10) to the photomicrograph data. For the basic 
fiber properties assumed, i.e.. I/d = 10, fiber eccentricity 

= 1.0, fiber modulus = 1.378 x 10 5 MN/m 2 (20xl0 6 PSI), fiber 

strength = 344.5 MN/m 2 (50KSI) , the figures indicate that the 
RSI material should have a modulus of elasticity between 1350 

2 ? 
and 1450 MN/m , and a strength between 2.5 and 2.7 MN/m , and 

a strain at the onset of instability of roughly 0.178%. 

Naturally, other assumptions regarding fiber modulus and strength 

would alter these results. A change in fiber SL/d value will alter 

the result, a small amount, similar to the changes seen in Figures 

21 and 22, but a change in fiber eccentricity would greatly alter 

the results as discussed in the following section. 

5.2.5 Effect of Fiber Eccentricity 

Initial fiber curvature affects were examined by 
considering fiber distributions with 0=0, distribution eccen- 
tricities of 0.0 and 0.3, and fiber length to diameter ratios 
of 10 and 15. The initial fiber curvature, fiber eccentricity, 
was varied from 0.0 to 3.0. The results are shown in Figures 
26 through 28. 

As can be seen in Figures 26 through 28, the fiber 
eccentricity produces the greatest overall effect on the strength 
characteristics of the material. The results show that decreas- 
ing the fiber eccentricity from 1.0 to near 0, one can obtain 
an increase in elastic modulus of 160%, a corresponding increase 
in strength of near 280%, along with 40% increase in strain to 
failure. Increasing the initial fiber eccentricity over 1.0 
results in decreases in modulus of elasticity and ultimate 
strength, but increases in strain to failure. The significant 
result here is shown in Figure 28 , in that the strain at first 
fiber break goes through a minimum value at slightly over an 
initial eccentricity of 0.50. 

A minimum strain value for a given fiber eccentricity 
indicates that if one desires to increase the strength of the 
material, processing which does not impart initial curvature 
of the fiber between intersections will produce the best results. 
Very small initial curvature will have a large effect on all 
properties. Large values of initial curvature will result in 
increasing strain to failure, but at the expense of both strength 
and modulus. Since a change in modulus will result in changing 
the thermal strains induced in the material, a trade-off study 



of thermal strain versus modulus and potential strain to failure 
versus modulus is required to determine the potential benefit 
of trying to alter the initial fiber curvature. 


6. DISCUSSION 


The previous sections have described the method of analysis 
and various parametric evaluations of RSI material which have 
been performed. The merit of the analytical model is apparent 
in that one can obtain, at a minimum, some "feel" for the type 
and magnitude of change in properties which one might expect from 
control of the various parameters during material manufacture. 
Thus, the analytical studies may be used as a basis for deter- 
mining the desired microstructural changes in RSI material to 
improve the mechanical properties of the material. 

Further confidence in the ability of the analytical model 
to predict RSI material properties can be obtained from compari- 
son of the predicted properties with those measured and reported 
on the material. A direct comparison is not easy however, since 
the fiber length to diameter ratio and initial curvature could 
not be easily defined from the photomicrographs. Consequently, 
a complete geometrical description of the material is not avail- 
able. In addition, the measured strength parameters exhibit con- 
siderable variations, further complicating a direct comparison 
between measured and predicted properties. It would appear that 
a reasonable comparison can be obtained if for the experimental 
range in modulus of elasticity, one chooses a range in model 
parameter (fiber eccentricity for example) to correspond to the 
same range of predicted modulus. Then, using the same model 
geometry ranges compare the expected range of ultimate tensile 
strength and strain to failure with the experimentally observed 
ranges . 

The above described approach has been used to generate 
Figure 29. The mechanical properties were those reported for 
the RSI-LDF. As can be seen, a value of fiber eccentricity 
ranging from 2.2 to 3.2 produces predicted elastic moduli in 
the measured range. Using the same range for fiber eccentri- 
city, the analytical model slightly over predicts the stress 
at failure but gives a reasonable estimate of the strain to 
failure. The analytical model is thus capable of predicting 
realistic mechanical properties for the RSI material. Veri- 
fication of the changes in properties with changes in micro- 
structure will naturally require extensive test data on samples 
with known microstructure. 

The results of the RSI material model studies indicate that 
small diameter fibers can provide equivalent structural capabil- 
ity to large diameter fiber structures , provided that the joint 
spacings are comparable. Structural properties have been shown 
to be a strong function of the orientation and eccentricity of 
the fiber distribution function, which indicates one possible 
means of tailoring material properties. In addition, initial 
fiber curvature has been shown to be a significant variable affec- 
ting properties, and hence, may provide a means of increasing 
strain to failure of future materials. 
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CONCLUSIONS AND RECOMMENDATIONS 


A suitable micromechanics model has been developed to assess 
the strength characteristics of rigidized fibrous material having 
a random oriented geometry in three dimensions. The basic assump- 
tions and the method of analysis are reported in section 
three of the present report. The merit of the model lies in its 
ability to deal with any arbitrary geometrical configurations 
and also any combination of mechanical and thermal loadings. The 
model can also deal with the spatial variability of other para- 
meters , like fiber eccentricity, length to diameter ratio, Young's 
modulus, strength of the fiber, etc. Results have been obtained 
for a large number of cases under different combinations of 
thermal and mechanical loadings. Some of these results when com- 
pared with experimental results obtained previously [6] for RSI 
material showed good agreement. 

It is naturally desirable to have more information on the 
experimental work done on different low density materials under 
different loading conditions. Comparison of these experimental 
works with the present analytical results will give a better under- 
standing of the accuracy and reliability of the present model. In 
this respect, we can^ say that based on the results obtained, a 
qualitative assessment of the merit of the present model has al- 
ready been established. 

The analytical model for RSI material has demonstrated 
sufficient accuracy in predicting properties. Consequently, future 
analytical studies should be based on exercising the present model 
with the aim of determining the material properties such as fiber 
distribution eccentricity and orientation, initial fiber curvature 
and material density which would be required to meet given sets 
of strength goals. Studies of this type should include goals such 
as combined loading strengths and combined mechanical and thermal 
stress capability. Thus, the analytical model can help guide 
future material development. 

Most of the results presented here have considered different 
cases of static loads and thermal loads. However, structural 
material is often subject to acoustic, vibratory or shock loads. 
Therefore, a study of their effects should be considered as a 
part of the future recommendation for research activity. 
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SYMBOLS LIST 


A - Cross sectional area of fiber 
a,b, c- Cell edge dimensions 
B - Transformation Matrix 

C* - Composite Stiffness Matrix 

d - Fiber diameter 

e - Strain components 

E - Axial fiber modulus 

E' - Effective fiber modulus 

l - Fiber length between intersections 
L - Total fiber length in cell 

n - Number of bays in cell 

v - Cell volume 

B^ - Orientation of fiber in given plane 

6 - Initial midspan fiber deflection 

6e - Variation in cell strain 

AE - Increase in cell strain energy 

- Orientation of vertical plane in cell 
a.. - Stress components 
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APPENDIX A - GOVERNING EQUATIONS 


A. Derivation of Stress/Strain Relation: 

The basic assumptions underlying the derivation of 
stress/strain relation of the space frame model have been dis- 
cussed in Section 2. With these assumptions, the strain com- 
ponent in the fiber direction of any truss is given by 


where 
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<fh is the angle of orientation of the truss with reference 

to the x axis, and 3-j_j is the angle of orientation of the 

fiber in any truss with the angle of orientation with 

reference to the x axis. Similarly, the stress in the fiber 
of any truss is given by 
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and in the case of compression in the fiber 
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At any stage of loading, the increase in strain energy of 
the model due to virtual displacement is given by 


AE = Ve 
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where Dij = £ B . E^j.B and the summation is carried out over 
all free lines in a J truss. Similarly, the additional external 
work done due to virtual displacement can be written as 

AW = V. a. 6e* (A-4 ) 
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Equating the total increase in strain energy to the additional 
external work for the equilibrium condition, we can write 


a = C * . e ; 
or a . . = c* 
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A. 2 Modification of Stiffness for Clamped Ends 

The effective value of Young's modulus used in the 
equation (A-2) represents the nonlinear axial stiffness for 
fiber links with simply supported ends. The necessary modi- 
fication for the clamped ends can be incorporated by applying 
moments of both ends such that the end slope becomes zero 
under axial compression or tension. In that case, the express- 
ions (A-2) for the nonlinear axial stiffness for any member are 
given by 
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for the case of compression, and by 
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for the case of tension. 

To prevent an unstable kinematic condition due to the 
relative rotation of bars of joints, the model must have 
at least one truss in the horizontal plane. 


APPENDIX B - RSI MODEL EVALUATIONS 

The basic formulation, failure criteria and. the method of 
analysis for the RSI model have been discussed in Section 3. 

The present section will deal with results obtained using the 
material model to study the effect of geometric and material 
variability on the strength of rigidized fibrous material unde- 
different loading conditions. The model essentially consists of 
a set of parallel trusses placed in five different planes as 
discussed in Section 3. The analysis was done initially con- 
sidering that the fiber between any two joints was free to ro- 
tate at the edge. However, subsequently the analysis has been 
modified to incorporate the reality that the fiber is clamped 
against rotation at both ends. Details of the results covering 
the strength of the model under tension, compression, shear and 
combined loadings are discussed below. 

B.l Comparison Between the Free and Clamped Edge: 

In Figure B-la, a comparison is drawn between the re- 
sults obtained for the free edge case and the clamped edge case 
The model considered consists of six trusses in each plane and 
is subject to pure tension in the vertical direction. Angle 
of orientation of the fibers and the variation of strength from 
truss to truss in any plane are shown in Figures B-lb and B-lc, 
respectively. As shown in Figure B-la, there is a considerable 
difference in results between free edge and clamped edge, 
and therefore, justifies the modification to incorporate an 



approximation to joint fixity. Results in Figure B-la also 
show a progressive failure as discussed under failure criteria 
in the previous section. 


B.2 Tension 

B.2.1 Study of the Distribution of Angle of- Orientation: 

Figures B-2 - B-6 refer to the different plots of the 
same model having six trusses in each of the five planes. The 
material and geometrical properties of the model are shown in 
Figure B-2. The model is subject to a uniform tension in the 
vertical direction and the ultimate strength of the model is 
determined for a constant volume fraction of 0.167 percent for 
each truss with the different distribution of the angle of 
orientation as shown in Figures B-5a to B-5h. The results obtained 
are shown in Table B-I. It is clearly observed that the different 
orientation has a significant influence on the ultimate strength 
of the model. Stress/strain plots shown in Figures B-2 to B-4 
refer to the distribution of the angle of orientation of the fiber 
as shown in Figures B-5a, b, and c, respectively. Curves marked 
(1) show the results obtained using variable volume content for 
different trusses as shown in Figure B-6, and curves marked (2) 
refer to the case when the volume content is constant. 


B.2. 2 Study of the Effect of the Strength Distribution: 

Table B-I I presents the results of the same model with 
the angle of orientation of the fiber given by Figure B-5b to 
study the effect of different strength distribution as shown in 
Figures B-7a to B-7e. Cases 2, 3 and 4 exhibited the progressive 
type of failure postulated by the authors in a previous section. 
Table B-I I clearly shows that the distribution of strength of 
the fiber can have a considerable effect on the ultimate strength 
of the RSI model. 


B.3 Compression 

B.3.1 Study of the Distribution of Angle of Orientation: 

Details of the space frame model for the study of the 
compressive strength of RSI material are shown in Figure B-8. 

It consists of four vertical planes oriented at angles of 0, 

+45, -45, 90 and one horizontal plane. Each of the vertical 
planes has six trusses with all parameters constant except 
the angle of orientation of the fibers. The values of the con- 
stant parameters are also shown in Figure B-8. The model is sub- 
ject to uniform compressive load in the vertical direction and 
the failure load is obtained for a constant volume content of 
the fibers equal to 5% with the different distribution of the 
angle of orientation as shown in Figures B-5a to B-5h. The 



results are given in Table B-III. These results when compared 
with the ultimate loads obtained for the case of tension are 
presented in Table B-I show a little difference in failure 
loads. This difference may, however, be significant at higher 
aspect ratios between fiber joints. One important thing to be 
noticed when comparing the results is that the effect of the 
distribution type of the angle of orientation of the fibers on 
failure loads is similar. 

B.3.2 Study of the Effect of the Strength Distribution: 

To observe how the variability of the strength dis- 
tribution affects the failure strength of RSI model, the ex- 
ample under Case (b) in Table B-III with the distribution type 
of the angle of orientation given in Figure B-5b were studied 
using different strength distributions as shown in Figures B-7a 
to B-7e. In all cases, trusses with similar fiber orientation 
in each plane have the same strength. Results obtained are 
tabulated in Table B-IV. It is clearly observed from Table B-IV 
that the distribution of strength has a significant effect on the 
ultimate strengths of the model, the nature of effect being the 
same as in the case of tension. Cases 2, 3 and 4 of Table B-IV 
exhibited the progressive type of failure as postulated by the 
authors under the combined effect of variability of the angle 
of orientation and strength distribution. 

B.4 Shear 

B.4.1 Study of the Distribution of Angle of Orientation: 

The parametric details of the space frame model used 
for the study of the shear strength of RSI material are shown 
in Figure B-8. It is the same model that is used for the study 
of failure loads under tension and compression. The model is 
subject to pure shear load in the 1-3 plane and the failure load 
is obtained for a constant volume content of the fibers equal 
to 5% with the different distribution of the angle of orienta- 
tion as shown in Figures B-5a to B-5h. Results obtained are 
given in Table B-V. These results when compared with the _ 
corresponding values obtained under pure tension show a marked 
difference. Whereas in the cases (a), (c) to (f) and (h) , 

failure loads obtained are lower, cases (b) and’ (g) show higher 
values for ultimate loads. The reason may be attributed to 
the different nature of effect of distribution angle on the 
ultimate strength of the material. 

B.4. 2 Study of the Effect of Strength Distribution: 

In the examples discussed above, strength of the fiber 
is assumed constant from truss to truss. However, in actual 
material, a sample of fibers shows a marked difference in their 
strength from fiber to fiber. To study this effect on the ulti- 
mate strength__of the model. Case (b).-in Table- B-V -is subject to 
different kinds of strength distribution as shown in Figures 
B-7a to B-7e, keeping the angle of orientation of the fibers the 
same. It should be noted that trusses in different planes with 
the same angle of orientation have, the same strength values; re- 
sults are given in Table B-VI. Also shown are the results obtained 
by reversing the order of maximum allowable stresses in the trusses 
in each plane. That is, since the initial model assigned a speci- 
fic fiber strength to each fiber, a different material model can 



be generated by simply reversing the strength assignments. The re- 
sult could be, for example, an initial case where the strongest 
fibers are oriented towards the applied load direction and a sub- 
sequent case where by simply changing strength assignments, the 
weakest fiber could be oriented towards the loading direction. It 
should be observed that changing the order of maximum allowable 
stresses from truss to truss, results for the ultimate stress of 
the RSI material can be significantly different. Because of these 
differences, one should model strength variability by creating N 
parallel trusses where each truss has its own strength. Thus, if 
M fiber angles are defined, one would create a model with N * M 
fiber pairs per plane. 

B.5 Combined Shear and Tension 

B.5.1 Study of the Distribution of Angle of Orientation: 

Details of the space frame model for the study of the 
strength of RSI material under combined tension and shear are 
shown in Figure B-8. It consists of four vertical planes oriented 
at angles of 0, +45, -45, 90 and one horizontal plane. Each of 
the vertical planes has six trusses with all parameters constant 
except the angle of orientation of the fibers. The values of 
the constant parameters are also shown in Figure B-8. The model 
is subject to combined uniform tension and shear in the vertical 
direction and the failure load is obtained for a constant volume 
content of the fibers equal to 5% with the different distribution 
of the angle of orientation as shown in Figures B-5a to B-5h. The 
results are given in Table B-VII . These results, when compared 
with the results obtained for the case of tension only and shown 
in Table B-I, show a significant drop in the ultimate stress due 
to the presence of shear. A similar comparison results when 
Table B-VII is compared with results obtained for the case of 
shear and shown in Table B-VI. 

B.5. 2 Study of the Effect of Strength Distribution: 

In the examples discussed above, strength of the fiber 
is assumed constant from truss to truss. However, in actual mate- 
rial, a sample of fibers shows a marked difference in their strength 
from fiber to fiber. To study this effect on the ultimate strength 
of the model. Case (b) in Table B-VII is subject to different kinds 
of strength distribution as shown in Figures B-7a to B-7e, keeping 
the angle of orientation of the fibers the same. It should be 
noted that trusses in different planes with the same angle of 
orientation have the same strength value. Results are given in 
Table B-VIII. Results clearly show that using the variable strength 
distribution, the ultimate stress of the material can be signifi- 
cantly reduced. 


APPENDIX C - RSI MATERIAL REPRESENTATION 


C.l Discrete Representation of a Continuous Distribution Function: 

The space frame analysis program requires representing con- 
tinuous distribution functions with a discrete number of points. 
Spatial orientation of fibers, fiber diameters, length between 


25 



fiber intersections, fiber eccentricities and fiber strengths 
all possess continuous distribution functions. Analytical 
representation of RSI material requires identifying a deter- 
ministic set of fiber properties which can be modeled with a 
finite number of trusses with fixed properties. The following 
section describes how a continuous distribution can be repre- 
sented by a discrete distribution of finite number of samples. 

In particular, the spatial distribution of fibers in any given 
plane is considered. However, once the continuous distribution 
function for any of the variables is defined, the discrete repre- 
sentation of the variable can be determined in an analogous 
fashion. 

It is first assumed that the continuous distribution func- 
tion of a variable is known, in this case the angular orienta- 
tion of RSI fiber in a given plane. It has been assumed that 
the distribution can be represented as a first approximation 
by an elliptical distribution f(9) such that 

f(0) = i (1+e cos 26) (C-l) 

where e is a measure of the eccentricity of the ellipse and 
6 is understood to be the angle of a fiber as measured from 
the mean angle or preferred fiber orientation. In practice the 
mean angle, 0, and eccentricity can be determined by using 
scanning lines, on photomicrographs of RSI material. Namely, 
the number of fiber intersections with the scanning line, 
versus orientation of the scanning line will provide an experi- 
mental distribution of f(9). A large number of scanning lines 
located at various areas of RSI material will be employed to 
give a good representation of f(0). From the experimental 
representation of f(0), the mean value of 9 can be found and 
the eccentricity can be determined in a least squares sense, 
the results of which are described in the following section. 

The discrete representation of f ( © ) is given by 
N/2 

g(<p) = I ^ 4 >i (C-2 ) 

-N/2 


where N is the total number of points and = implying 

that <f>^ is normalized with respect to the mean such that the 
mean values of f(9) and g (<p) are zero. Specific values of 
are determined as follows: 

into N regions such that 

n = 0 , 1 , . . . , N (C— 3) 


f(0) is divided 
6 

/ 0 n f (9 ) d9 = 


26 



and one value of <{k is located within each resulting interval. 

Each <b. could be placed anywhere within the interval 0 . to 

1 n- 1 

0 , and the general solution requires satisfying the first N+l 
moments of f (9) in order to determine each (Jk location. How- 
ever, a reasonable approximation of f(0) can be made by placing 
all values of <j>^ in the same relative position within 8 n _-^ and 

0 n and minimizing the difference between the second moments of 

f(9) and g (<}>). One thus obtains the following expression to 
be solved for a, the position of the discrete variable <}>^ 

between 0 and 0 
n-1 n 


N/2 

Z 

1 


2 

N 


{0 


N-1 


+ a 


( 0 . 


9 „-i> > 2 - s ■ 0 


(C— 4 ) 


where S = i (?■ - e) is the variance of f(9) 

Z o 

0 . ■* 9 defines the equal area intervals of f(9) 

n-1 n ^ 

N is the number of fibers = j the number of dis- 
crete trusses (each truss has +9 fibers) 


An example is tabulated in Table c-1 giving the desired 
positions of 4^ for various values of N. A graphical repre- 
sentation of two continuous distributions for e=0.3 and 0.95 

along with the discrete representations of the distributions 
are given in Figures C-1 and C-2. 

The above approach is valid for a standard elliptical 
distribution. However, when one considers how a preferred 
mean value of fiber orientation may occur during manufacture, 
a slight modification to the elliptical distribution appears 
appropriate. Assuming that the manufacturing process starts 

with a slurry of completely randomly oriented fibers, and 

then pressure is applied to compress the slurry, it is 
reasonable to assume that the fiber distribution will be 
symmetric about both the in-plane and through-the-thickness 
directions. This can be accomplished by defining f(9) such 
that 

f ( 0 ) = const. (1 + e cos 2 (0 - j0) ) (C-5) 

where 

j = 1 if 0 < 6 £ tt/2 

it < 9 <_ 3ir/2 

j = -1 if tt/2 < 9 <_ tt 
3tt/2 < 9 < 2it 
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and results in the distribution shown in Figure C-3. The 
function, f(0), must be normalized such that 

/ 1T f (0 ) d0 = 1.0 (C-6) 

which results in 

f( 9 ) = 1 + e cos 2 (0- j0 ) (c _ 7) 

■n + 2e sin 20 

Note that, when normalized, the constant for f(0) becomes 
a function of the orientation and eccentricity of the dis- 
tributor. Although the process is slightly more involved 
algebraically, equal area intervals and discrete locations 
for placement of discrete points can be determined in the 
same manner as previously outlined. In addition, data points 
taken from RSI photomicrographs can be fit to equation c-7 
in the least squares sense to define 0 and e ; the solution 
involving solution of a pair of non-linear equations to deter- 
mine 0 and e. 


C . 2 Least Squares Curve Fit to Distribution Functions 

Given data from the photomicrographs of RSI material which 
can be tabulated as number of fiber intersections, f^, with 

scanning line versus angular orientation, 0^, of scanning line, 

one can determine the best fit values of e and 0 in the follow- 
ing manner. 

Define the sum of the squared residuals as 

V = ? {f. - 1 + e cos 2 (9_- 8_j } 2 (c _ 

1 it + 2e sin 2 0 

desired to minimize V by proper selection of 
require 


(C-9) 


Taking the partial derivatives, and after some manipulation, 
one can solve each equation explicitly for e in terms of 0. 
as shown below: 


Since it is 
e and 0 , we 


3V 

3e 


= 0 


9 V 
30 


= 0 



From 


tt— = 0 one obtains 
9e 

A 2 cos 20 - A^ sin 2 0 

0 o — _ — 

A 0 sin 20 + A. cos 29 + A,- sin 20 cos 20 
3 4 5 

where 

A x = 2tt E - 2N - it £ (irf i ~l) sin 2© i 
A 2 = - it E (irf^-1) cos 20^ 

A, = 4 E f. + tt * E sin^ 20. - 2E (Tif.+l) sin 20. 

3 x x x x 

2 

A^ = tt E cos 2 0^ 

A c = 2E sin 20. cos 20. - 2 E (uf.+l) cos 20. 

5 x x x x 

and where E represents summation over the number of experi- 
mental data pairs, N. 

A similar, but more lengthy, expression is obtained 

from 


80 

Solution can be obtained by simply substituting various values 
of 0 and plotting the resultant values of e until an inter- 
section of the two curves is obtained, noting the restrictions 
that the derivation requires e <_ 1.0, and 0 £ 0 £ tt/2. 


C.3 Results Obtained From RSI Photomicrographs 

The photomicrographs supplied to Materials Sciences Cor- 
poration by NASA/Langley have been studied and data taken to 
provide fiber distribution information. The scanning line 
technique was used as follows. A line 600 micrometers in length 
was taken as a fixed length. Angular orientations from -90 
to +90 degress, in 15 degree increments were marked on the 
photographs. For each angular orientation, the scanning line 
was placed at twenty different locations on the photograph 
and the number of fibers which intersected the scanning line 
were counted. Thus, for each of twelve angular orientations, 
twenty readings of number of intersections were obtained in 
an effort to obtain representative numbers from each photo- 
graph. Fibers which could be identified through the stereo 
photos as being significantly out of the plane of the photo- 
graph were not included when totaling the number of fiber in- 
tersections. Since the angle which yields the most inter- 
sections implies that the primary fiber orientation is at 
ninety degrees to that angle, all angles were shifted ninety 
degrees to reflect the actual fiber orientation. 



The resulting tables of data were then used to determine 
the best fit values of e and 9 for equation C-l. The results 
obtained from equation c-l along with considering the fabri- 
cation process prompted fitting the data to equation c-7 since 
it was felt that this distribution function would be more 
realistic. 



RSI MODEL RESULTS FOR TENSION LOADING 
EFFECT OF ANGLE OF ORIENTATION OF FIBERS 


TABLE B-I 

Ultimate Stress of 

Distribution Type of the the RSI Material 

Angle of Orientation of 2 

the Fibers MN/m (PSI) 


Fig . 

B-5 (a) 

1.405 

(204) 

II 

" (b) 

1.275 

(185) 

II 

" (c) 

1.536 

(223) 

II 

" (d) 

1.902 

(276) 

it 

" (e) 

1.006 

(146) 

II 

" (f) 

2.150 

(312) 

M 

(g) 

1.226 

(178) 

II 

" (h) 

1.330 

(193) 


RSI MODEL RESULTS FOR TENSION LOADING 
EFFECT OF FIBER STRENGTH DISTRIBUTION 


TABLE B- II 


Distribution Type 
of the Angle of 
Orientation of the Fiber 

Strength 
Type of 

Distribution 
the Fiber 

Ultimate Stres 
of RSI Model 

MN/m 2 (PSI) 

Fig. B-5 

(b) 

Fig . 

B-7 

(a) 

1.275 

(185) 

II II 

(b) 

II 

II 

(b) 

.661 

( 96) 

II II 

(b) 

II 

II 

(c) 

.696 

(101) 

II II 

(b) 

II 

•1 

<d) 

.661 

( 96) 

tl II 

(b) 

II 

II 

(e) 

. 785 

(114) 
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RSI MODEL RESULTS FOR COMPRESSION LOADING 
EFFECT OF ANGLE OF ORIENTATION OF FIBERS 


TABLE B-III 


Distribution Type Of The 
Angle Of Orientation 
Of The Fibers 


Ultimate Stress of 
the RSI Material 

MN/m 2 (PSI) 


Fig. B-5 


(a) 

(b) 

(c) 

(d) 

(e) 

(f) 

(g) 

(h) 


-1.350 

(-196) 

-1.233 

(-179) 

-1.481 

(-215) 

-1.826 

(-265) 

- .971 

(-141) 

-2.067 

(-300) 

-1.213 

(-176) 

-1.281 

(-186) 


RSI MODEL RESULTS FOR COMPRESSION LOADING 
EFFECT OF FIBER STRENGTH DISTRIBUTION 


TABLE B-IV 


Distribution Type Of The 
Angle of Orientation 
Of The Fibers 

Strength 
Distribution 
Type of Fiber 

Ultimate 
the RSI 

MN/m 2 

Stress 

Materi; 

(PSI) 

Fig. B-5 (b) 

Fig. B-7 

(a) 

-1.233 

(-179) 

" " (b) 

11 11 

(b) 

- .648 

(- 94) 

" " (b) 

II II 

(c) 

- .668 

(- 97) 

" " (b) 

II II 

(d) 

- .586 

(- 85) 

" " (b) 

It II 

(e) 

- .772 

(-112) 
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RSI MODEL RESULTS FOR SHEAR LOADING 
EFFECT OF ANGLE OF ORIENTATION OF FIBERS 


TABLE B-V 


Distribution Type Of 
The Angle Of Orientation 
Of The Fibers 


Ultimate Stress of 
the RSI Material 

MN/m 2 (PS I) 


Fig . 

B-5 (a) 

1.254 

(182) 

II 

" (b) 

1. 357 

(197) 

II 

" (c) 

1.302 

(189) 

II 

" (d) 

. 572 

( 83) 

II 

" (e) 

.854 

(124) 

II 

" (f) 

.675 

( 98) 

II 

" (g) 

1.709 

(248) 

II 

" (h) 

.999 

(145) 


RSI MODEL RESULTS FOR SHEAR LOADING 
EFFECT OF FIBER STRENGTH DISTRIBUTION 


TABLE B-VI 


Distribution 
Type of the 
Angle of 
Orientation 
of the Fibers 


Strength 
Distribution 
Type of 
Fiber 


Ultimate Stress of 
the RSI Material 


MN/m 2 

Results For Initial 
Fiber Strength 
Assignments 


(PSI) r 

Results For Revising 
Order of Fiber Strength 
Assignments 


(b) 

Fig.B-7 

(a) 

1. 357 

(197) 

1.357 

(197) 

(b) 

II II 

(b) 

1.158 

(168) 

.937 

(136) 

(b) 

II II 

(c) 

1.089 

(158) 

.923 

(134) 

(b) 

II II 

(d) 

1.075 

(156) 

.902 

(131) 

(b) 

II II 

(e) 

1.309 

(190) 

1.054 

(153) 
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RSI MODEL RESULTS FOR COMBINED 
TENSION AND SHEAR LOADING 
EFFECT OF ANGLE OF ORIENTATION OF FIBERS 


TABLE B-VII 



Distribution Type 
Of The Angle Of 
Orientation Of Fibers 

Ultimate 
the RSI 

MN/m 2 

Stress of 
Material 

(PSI) 

Fig 

. B-5 (a) 

.827 

(120) 

If 

II 

(b) 

.820 

(119) 

II 

„ (c) 

.84 7 

(123) 

II 

- (d) 

.482 

( 70) 

If 

" (e) 

.537 

( 78) 

II 

" (f) 

.558 

( 81) 

II 

" (g) 

1.206 

(175) 

IT 

" (h) 

1.24 7 

(181) 


RSI MODEL RESULTS FOR COMBINED 
TENSION AND SHEAR LOADING 
EFFECT OF FIBER STRENGTH DISTRIBUTION 


TABLE B-VII I 


Distribution Type 
Of The Angle Of 
Orientation Of Fibers 


Strength 

Distribution Type 
Of The Fiber 


Ultimate Stress of 
the RSI Material 

MN/m 2 (PSI) 


Fig . 

B-5 (b) 

Fig. 

B-7 (a) 

. 820 

(119) 

II 

" (b) 

H 

" (b) 

.441 

( 64) 

It 

II 

it 

II 




(b) 


(c) 

.427 

( 62) 

II 

" (b) 

it 

•• (d) 

.400 

( 58) 

II 

" (b) 

it 

" (e) 

.530 

( 77) 
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RSI MODEL RESULTS FOR THERMAL LOADING ON MODEL IN FIG. 13, 
GIVEN INITIAL THERMAL STRAINS AND APPLIED TENSION LOADING 

TABLE B-IX 


Distribution Type 
Of The Angle Of 
Orientation Of Fibers 


Initial Strains e°^ And £° 

Applied In The 1 And 2 
Directions 


Ultimate Stress of 
RSI Material In Ten- 

sion MN/m^ (PSI) 


Fig . 

B-5 (b) 

0 

1.2 75 

(185) 

II 

" (b) 

0.000424 

1.254 

(182) 

II 

" (b) 

0.000784 

1.240 

(180) 

II 

" (b) 

0.001144 

1.226 

(178) 

II 

" (b) 

0.001504 

1.213 

(176) 

II 

" (b) 

0.001864 

1.199 

(174) 


RSI MODEL RESULTS FOR THERMAL LOADING ON MODEL IN FIG. 13, 

GIVEN INITIAL THERMAL STRAINS AND COMBINED TENSION AND SHEAR LOADING APPLIED 

TABLE B-X 


Distribution Type 
Of The Angle Of 
Orientation Of Fibers 


Initial Strains And e 

Applied In The 1 And 2 
Directions 


Ultimate Stress of 
RSI Material In Ten- 

sion MN/m 2 (PSI) 


Fig . 

B-5 (b) 


II 

" (b) 

0 

II 

" (b) 

0 

II 

” (b) 

0 

II 

" (b) 

0 

II 

" (b) 

0 


0 

.820 

(119) 

000424 

.806 

(117) 

000784 

.785 

(114) 

001144 

.772 

(112) 

001504 

.751 

(109) 

001864 

.730 

(106) 
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TABLE C-I 


Values of for Representation of f(9) at e <= 0.95 
1 . * 


2 

4 

6 

8 

33.96 

8.52 

5.54 

4.19 


47.26 

22.44 

16.32 



54.08 

30.37 

58.36 


.3773 .3467 .3502 .3569 


10 

12 

16 

20 

3.40 

2.88 

2.22 

1.93 

12.94 

10.75 

8.08 

6.60 

23.30 

19.07 

14.10 

11. 36 

35.72 

28.34 

20.47 

16.28 

61.38 

39.70 

27.41 

21.45 


63.65 

35.35 

27.01 



45.34 

33.18 



66.93 

40.35 




49.56 




70.11 

.3639 

.3703 

.3819 

.4156 




A) Small Diameter Fiber Material - RSI-SDF 
Magnification: 200 X 



B) Large Diameter Fiber Material - RSI-LDF 
Magnification: 100 X 


Figure 1- Typical RSI Photomicrographs 
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LAYOUT OF VERTICAL TRUSSES 




0 


TYPICAL VERTICAL TRUSS 
n = NUMBER OF "BAYS" 

= NUMBER OF SEGMENTS 
PER FIBER 
L= FIBER LENGTH 
- LAYOUT OF HORIZONTAL TRUSS 


Figures - A typical space model. 
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.07 


a j = 0, + 45, -45, 90 
Bj 3 60° IN ALL CASES 
i/d =10, Ejj =20x|0 6 
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MN/m 


TRUSS E 
NO MN/m 2 


F 

MN/m 2 


®= 0.689 xIO 5 = 103.3 

©= 1.034xl0 5 = 206.7 

®= 1.722 xIO 5 = 206.7 



Figure 5 - Stress-strain behavior of the model under tension and 
compression . 
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COATED SURFACE 



Figure 5 - Position and orientation of SEM specimens in contractor 
supplied RSI tiles. 



'-f(0)»i I +.188 COS 2(0-175°) 
If ' 


9 = 17.5° 
e =.188 


Figure 7 - Fiber distribution data for rsi -LDP - X-Z Plane 








9 = -87.85° 
e= 5.81 x 10" 5 

i.e DENSITY FUNCTION 
IS CIRCULAR 


Figure 10 - Fiber distribution data for RSI-SDF - X-Y Plane 



, <nx l+e COS ( 9 - 9 ) 

f (a) 7* 

1T+2e SIN 2 0 


I 


where 9 = 18° 

e =0.070 


Figure 


11 - 


Fiber distribution 
distribution for 


data and best fit elliptical 
RSI -LDF - X-Z Plane 



*Y 



f(0) = 


Ue COS 2(0-0) 
IT- 2e SIN 2 0 


where 


0=23.4° 


e = 0.148 


Figure 


2_2 ~ Fiber distribution data for best fit elliptical 
distribution for RSd-LDF - Y-Z Plane 
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AN 
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Figure 14 - Fiber Diameter Variation 

a) RSI-SDF X-Y Plane 

b) RSI-LDF - X-Y Plane 


PLANE 



Figure ^5 - Fiber Diameter Variation 

a) -RSI-LDF - Y-Z Plane 

b) RSI-LDF - x-Z Plane 
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PLANES 0,±45,90 
LOADING cr 33 

FIBER DISTRIB. e=0.0 
FIBER ECC.= 1.00 
L/D® 5,10, 15 
FIBER UTS =344.5 MN/m 2 
V f =O.IO 


a : 

UJ 

mcvi 


L/D = 



Figure 15 - Effect of the number of fiber pairs/plane on 
analysis results. 
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PLANES 0, ±45,90 
LOADING <r 33 

FIBER DIST. e =0.2 
FIBER ECC = 1.00 
L/D = 10 

FIBER UTS =344.5 MN/m 2 
V f = 0.10 


oc 



=5xl0 6 m =3.l8xl0" 3 m fl£,L/D) 
CASES 


Figure 20 


E-ffect of varing fiber diameter while holding 
fiber L/D constant. 
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MODULUS @ I FIBER STRAIN @ I st FIBER 

BREAK - MN/m 2 BREAK % 


PLANES 0,i45 t 90 
LOADING o- 33 

FIBER ECC.= 1.00 
L/D = 5, 10, 15 
FIBER UTS * 344.5 MN/m 2 
# FIBER PAIRS/PLANE =8 
V f = 0.10 

L/D =15 

^ • — — — ■ mmm m mm • 

. 10 

■ • ■' " ' ■ • 

“ “ • - ■ ■ • 

5 

i i l I I i i j i 



ECCENTRICITY OF FIBER DISTRIBUTION 



Figure 21 - Effect of eccentricity of fiber distribution 
on properties, 9 = 0.0 


0. 



UTS 33 STRESS AT I st FIBER BREAK 
MN/m^ 


. PLANES 0,* 45,90 
LOADING o - 33 
FIBER ECC.= 1.00 



ECCENTRICITY OF FIBER DISTRIBUTION 


Figure 22 ~ Ef feet _of -eccentricity of -fiber distribution 
on strength, 9 =0.0 



MODULUS AT I 5 *' FIBER BREAK E,, MN/m 


PLANES -0,^45, 90 
LOADING a 33 

L/D* 10 
FIBER ECC = 1.0 
FIBER UTS =344.5 MN/m 2 
V f =0.10 


= 0.3 


N0TE : "O" INDICATES LOCATION 
OF "BEST FIT" DATA TO 
RSI PHOTOMICROGRAPHS 
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ORIENTATION OF MEAN FIBER ANGLE 
(DEGREES) 


Figure 23 -■ Effect of orientation of fiber distribution on 
modulus of elasticity. 




ORIENTATION OF MEAN FIBER ANGLE 
(DEGREES) 

Figure 24 - Effect of orientation of fiber distribution on 
strength. 
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STRAIN AT I 8 ' FIBER BREAK 


PLANES - 0, ±45,90 



ORIENTATION OF MEAN FIBER ANGLE 
(DEGREES) 


Figure 25 ~ Effect of orientation of fiber distribution 
on strain at first fiber break. 
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FIBER ECCENTRICITY 

Figure 26 - Effect of fiber eccentricity on modulus of elasticity. 






PLANES - 0,* 45, 90 
LOADING o- 33 

FIBER DISTRIBUTION e =0.0,03 0 = 0.0 
L/D = 10, 15 

FIBER UTS = 344.5 MN/m 2 



0.5 


FIBER ECCENTRICITY 


Figure 28 “ Effect of fiber eccentricity on strain at first 
fiber break. 
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ECCENTRICITY FROM 2.2 TO 3.2 









Figure B-la - Stress-strain behavior of RSI material considering 
strength variability (volume fraction = 5%). 






STRESS, 0-a MN/m 


MATERIAL AND GEOMETRICAL PROPERTIES 
a,* =0,+45, -45,90 



STRAIN e 33 


Figure B-2 - Stress-strain behavior of RSI model under tension. 
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STRESS, MN/m^ 






ANGLE OF ORIENTATION OF FIBER (ft) ANGLE OF ORIENTATION OF FIBER (ft) ANGLE OF ORIENTATION OF FIBER (ft) ANGLE OF ORIENTATION OF FIBER (ft) 



2 3 4 5 6 

TRUSS NUMBER 



TRUSS NUMBER 




TRUSS NUMBER 

Figure B-5 - Variability of angle of orientation of the fiber 
different trusses in any plane. 
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VOLUME FRACTION 


Figure B- 



TRUSS NUMBER 


- Variability of volume fraction of different trusses 
in a plane. 






MATERIAL AND GEOMETRICAL 


PROPERTIES' 


L/D = 10* 
e = 1.0; 


n = 4.0* 

L = 1.0* 

E f * J.378xl0 5 MN/m 
F =344.5 MN/m 2 


2 


SAME 
* FOR ALL 
TRUSSES 


j 


= 0, +45,-45,90 
LAYOUT OF VERTICAL TRUSS. 
SIX TRUSSES IN EACH PLANE 


Figure B-8 - RSI model used for compression, shear and 
combined tension/shear loading. 








_ I + e COS 2(0- j 9 ) 

“ 7T+ 2e SIN 2 £ 

j = + I if CK0<ir/2 
7T< 0 < 3 7T/2 
j = - I if 7r/2 < 0 < 7T 
<9 <Ztr 


Figure C-3 - Fiber distribution for RSI material 


NASA-Langley, 1974 
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